Report for Model 3

— for both jointly and separately

2024-08-10

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.00
## LDevsum       1.00       1.01
## dh0           1.10       1.26
## dh1           1.16       1.41
## dh2           1.05       1.13
## dh3           1.01       1.02
## dl0           1.04       1.06
## dl1           1.12       1.32
## 
## Multivariate psrf
## 
## 1.11

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1146.14620 22551.3845
DIC3 1112.94318 22642.5526
PWAIC 23.35143 206.4733
WAIC 1124.71772 22656.2759

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.03
## LDevsum       1.00       1.00
## dh0           7.43      13.83
## dh1          16.90      30.91
## dh2          14.34      32.04
## dh3           8.50      16.56
## dl0           1.25       1.63
## dl1           1.10       1.27
## 
## Multivariate psrf
## 
## 20.1

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1387.28489 33774.328
DIC3 1249.46624 33888.911
PWAIC 92.30564 167.363
WAIC 1319.99773 33893.412